3750-3767 Nucleic Acids Research, 2014, Vol. 42, No. 6 
doi:10.1093/}mr/gktl359 



Published online 13 January 2014 



The light-induced transcriptome of the zebrafish 
pineal gland reveals complex regulation of the 
circadian clockwork by light 

Zohar Ben-Moshe^'^, Shahar Alon^'^'^, Philipp Mracek^, Lior Faigenbloom^, Adi Tovin^, 
Gad D. Vatine\ Eli Eisenberg^ '*, Nicholas S. Foulkes^ and Yoav Gothilf^'^ * 

^George S. Wise Faculty of Life Sciences, Department of Neurobiology, Tel-Aviv University, Tel-Aviv 69978, 
Israel, ^Sagol School of Neuroscience, Tel-Aviv University, Tel-Aviv 69978, Israel, ^Institute of Toxicology and 
Genetics, Karlsrulne Institute of Teclnnology, Eggenstein-Leopoldshafen 76344, Germany and "^Raymond and 
Beverly Sacl<ler Sclnool of Physics and Astronomy, Tel-Aviv University, Tel-Aviv 69978, Israel 

Received November 4, 2013; Revised December 6, 2013; Accepted December 9, 2013 



ABSTRACT 

Light constitutes a primary signal wliereby en- 
dogenous circadian clocks are synchronized 
('entrained') with the day/night cycle. The molecular 
mechanisms underlying this vital process are known 
to require gene activation, yet are incompletely 
understood. Here, the light-induced transcriptome 
in the zebrafish central clock organ, the pineal 
gland, was characterized by messenger RNA 
(mRNA) sequencing (mRNA-seq) and microarray 
analyses, resulting in the identification of multiple 
light-induced mRNAs. Interestingly, a considerable 
portion of the molecular clock (14 genes) is light- 
induced in the pineal gland. Four of these genes, 
encoding the transcription factors dec1, reverbbi, 
e4bp4-5 and e4bp4-6, differentially affected clock- 
and light-regulated promoter activation, suggesting 
that light-input is conveyed to the core clock 
machinery via diverse mechanisms. Moreover, we 
show that ded, as well as the core clock gene 
per2, is essential for light-entrainment of rhythmic 
locomotor activity in zebrafish larvae. Additionally, 
we used microRNA (miRNA) sequencing (miR-seq) 
and identified pineal-enhanced and light-induced 
miRNAs. One such miRNA, miR-183, is shown 
to downregulate e4bp4-6 mRNA through a 3'UTR 
target site, and importantly, to regulate the 
rhythmic mRNA levels of aanat2, the key enzyme 
in melatonin synthesis. Together, this genome- 
wide approach and functional characterization of 
light-induced factors indicate a multi-level regula- 
tion of the circadian clockwork by light. 



INTRODUCTION 

Intrinsic circadian oscillators driving daily rhythms of 
physiology and behavior evolved as an adaptation to the 
24-h solar cycle (1). At the molecular level, the vertebrate 
circadian clock is based on transcriptional and transla- 
tional auto-regulatory loops that generate oscillations of 
clock genes and their protein products (2). In the primary 
feedback loop are the positive components, CLOCK and 
BMAL, that heterodimerize and activate transcription of 
target genes, including period (per) and cryptochrome {cry), 
through E-box enhancer elements. Negative feedback is 
accompHshed by PER:CRY heterodimers that re-enter 
the nucleus and repress their own transcription by inhibit- 
ing the CLOCK:BMAL complex (3). The robustness and 
accuracy of this core loop is achieved by stabilizing acces- 
sory feedback loops and post-translational regulation 
(3,4). Nevertheless, the intrinsic period length of the 
circadian clock deviates from 24 h and therefore requires 
daily synchronization by environmental time-cues, among 
which hght is the most dominant (5). Changes in gene 
expression within central clock structures is thought to 
be one of the processes required for Ught-entrainment, 
but what exactly these changes are and how they lead to 
hght-entrainment is yet to be defined. 

The zebrafish serves as an attractive vertebrate model 
for exploring the central circadian clock and its entrain- 
ment by light (6). As in other non-mammalian vertebrates, 
the zebrafish pineal gland incorporates all the elements 
required for light-entrainment and circadian rhythm gen- 
eration and is considered to serve as a master clock organ; 
it is photoreceptive and contains an intrinsic circadian 
oscillator that controls the rhythmic production of the 
hormone melatonin (7). Furthermore, most zebrafish 
tissues and even cell lines possess peripheral clocks that 
can be entrained by direct light exposure (8). This differs 
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from the case in mammals, in which the retina is the 
source of photic input to the clock, and the peripheral 
clocks are entrained indirectly by the master oscillator 
that is located in the suprachiasmatic nucleus (SCN) of 
the hypothalamus (9). Accordingly, the zebrafish pineal 
gland and zebrafish cell hnes constitute a suitable 
platform for studying the molecular mechanisms of the 
circadian clock and its entrainment by light. 

Light-induced gene expression constitutes a key element in 
the relay of photic information to the core clock machinery. 
In zebrafish, light exposure triggers the transcription of per2 
and cry la in most tissues and in cell lines, a process which is 
thought to underlie the entrainment of the circadian clock 
(10,11). Systematic promoter analyses of these clock genes 
identified D-box as a key cAv-acting element in the regulation 
of Ught-driven gene expression (12,13). Transcriptome 
analyses in zebrafish embryos (14), larvae, heart and cell 
cultures (15) demonstrated light regulation of various 
cellular processes, such as transcriptional control and 
DNA-repair, using microarray technology. These studies 
have contributed to the understanding of light-regulated 
gene expression in peripheral clock-containing tissues. 
Nevertheless, a comprehensive understanding of light- 
entrainment requires investigation of light-induced gene ex- 
pression in the central clock organ. Moreover, high-through- 
put messenger RNA (mRNA) sequencing (mRNA-seq) now 
offers improvements in the profiling of the light-induced tran- 
scriptome, mainly due to unlimited gene representation and 
generation of enhanced quantitative data (16). 

With the aim of further unraveling the mechanisms 
underlying light-entrainment of the vertebrate central cir- 
cadian clock, we used both mRNA-seq and microarray 
technologies to characterize the hght-induced transcrip- 
tome in the zebrafish pineal gland. Focusing on hght- 
induced genes that encode transcription factors {decl, 
reverhhl, e4bp4-5 and e4bp4-6), we demonstrated their 
effects on clock- and light-regulated promoter activation. 
Among these factors, we showed that decl, as well as the 
core clock gene per2, is important for the light-entrain- 
ment of rhythmic locomotor activity. Light-induced 
genes generally exhibit a transient increase in expression, 
pointing to the existence of mechanisms that govern 
niRNA stabihty, such as regulation by microRNA 
(niiRNA). Accordingly, we used miRNA sequencing 
(niiR-seq) to identify miRNAs whose expression is 
enhanced in the pineal gland and induced by light. We 
demonstrated that the pineal-enhanced and hght-induced 
miR-183 regulates e4hp4-6 mRNA levels through a target 
site in its 3'UTR. Moreover, we showed that miR-183 
contributes to the generation of rhythmic mRNA levels 
of the pineal clock-controlled gene, aaiiat2. Overall, the 
genomic approach and functional analyses indicate a 
complex regulation of the central circadian clock by light. 



MATERIALS AND METHODS 

mRNA-seq, miR-seq and microarray 
experimental procedures 

Adult (0.5-1.5-year-old) Tg[aanat2: enhanced green fluores- 
cent protein (EGFP)] zebrafish (17) were raised in a 
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Figure 1. Experimental design for transcriptome analysis. Horizontal 
bars represent the photic treatments preceding tissue sampling: 
groups of fish were exposed to a 1-h light pulse (white boxes) or 
kept in darkness (black boxes), and tissues were sampled immediately 
after the treatment at six circadian time points indicated by arrows. 



temperature-controlled re-circulation water system under 
12-h: 12-h light-dark (LD) cycles, and transferred to 
constant darkness (DD) at the end of the day before the 
experiment. Fish were exposed to a 1-h light pulse 
('daylight' fluorescent lamp, light intensity of 5W/m^) 
before sampling (light treatment) or kept under constant 
darkness for control (dark treatment). The fish were 
anesthetized in 1.5 mM Tricane (Sigma) and sacrificed by 
decapitation. Pineal glands and brains were removed and 
immediately processed within 3 min of dissection under an 
Olympus dissecting microscope SZX12 equipped with filters 
for excitation (460-490 nm) and emission (510-550 nm) of 
EGFP; the use of transgenic fish expressing EGFP in the 
pineal gland enabled its selective removal. The tissues were 
collected from light- and dark-treated fish at six time 
points with 4-h intervals throughout one daily cycle, corres- 
ponding to circadian times (CT) 2, 6, 10, 14, 18 and 22 
(Figure 1). 

For microarray analysis, 12 pineal glands were collected 
and pooled at each time point under both light and dark 
treatments. For mRNA-seq, pineal glands were collected 
at the six sampling times and pooled together to estabhsh 
triplicates of 15 pineal glands for every treatment. For 
miR-seq, 18 pineal glands were pooled from the six 
sampling times for both light and dark treatments, and 
one brain pool was composed of four brains that were 
coflected from hght- and dark-treated fish at CT2 and 
CT14. Brain samples were collected subsequent to the 
removal of pineal glands and therefore do not contain 
pineal tissues. Total RNA for mRNA and miRNA 
analyses was isolated using RNeasy Lipid Tissue Mini 
Kit (Qiagen) and miRNeasy Mini Kit (Qiagen), respect- 
ively. The sequencing data were deposited to the Sequence 
Read Archive, under accession SRP016132. The 
microarray data were deposited to the Gene Expression 
Omnibus, under accession GSE53288. 



3752 Nucleic Acids Research, 2014, Vol. 42, No. 6 



mRNA-seq and data analysis 

The lUumina TruSeq protocol was used to prepare libraries 
from RNA samples. Overall, six libraries (three light groups 
and three dark groups) were run on one lane of an lUumina 
HiSeq2000 machine using the multiplexing strategy of the 
TruSeq protocol (Institute of Applied Genomics, Italy). On 
average, ~6.7 miUion paired-end reads were obtained for 
each library. The reads were of 2 x 100 bp. TopHat (18) 
was used to align the reads against the zebrafish genome, 
keeping only uniquely aligned reads with up to two 
mismatches per read. Different pair-end reads with identical 
start positions in the genome in both pairs might represent 
polymerase chain reaction (PGR) duphcates (19); therefore, 
in these cases, only one representative pair-end read was 
included. Only reads that are aligned to the protein coding 
regions of known NGBI reference sequence (RefSeq) genes 
were used. A custom script written in Perl was used to parse 
the output of TopHat, given in Sequence Alignment/Map 
format (http://samtools.sourceforge.net/), and to convert it 
into raw number of reads aligned to each position in each 
ReBeq gene. The RefSeq genes information was obtained 
from the Table Browser of the UGSG genome browser (gen- 
ome. ucsc.edu/) using the zebrafish July 2010 (Zv9/danRer7) 
assembly. On average, ~10yo of the reads passed all filters 
and were mapped to coding regions of RefSeq genes. 

Statistically significant differences between the number 
of reads ahgned to each RefSeq gene (the expression 
profile) in the different tested conditions were identified 
as described (20). Briefly, the expression profiles were 
normalized using a variation of the trimmed mean of M- 
values normalization method (20,21). Subsequently, we 
searched for expression differences between the three 
dark samples and the three fight samples that cannot be 
explained by the expected Poisson noise with P < 0.05 and 
Bonferroni correction for multiple testing (20). 

miR-seq and data analysis 

miRNA capturing and hbrary construction for the pineal 
fight, pineal dark and brain samples were conducted using 
lUumina's TruSeq Small RNA Sample Prep Kit according 
to the manufacturer's protocol. The mature miRNA 
fibraries were sequenced with bar-codes on one lane of 
an lUumina HiSeq2000 instrument following the manufac- 
turer's protocol. The reads were filtered using lUumina's 
HiSeq2000 software. The total number of reads after 
filtering was ~6, ~10.4 and ~14.8 niilfion for the pineal 
fight, pineal dark and brain, respectively. Next, the reads 
were aligned against the zebrafish genome, as described 
(22), providing ~4, ~5.8 and ~4.3 milfion reads for the 
pineal fight, pineal dark and brain, respectively. 
Approximately 3.8, ~5.4 and ~3.5 million reads were 
aligned to known miRNAs [miRBase, release 17; (23)] 
for the pineal light, pineal dark and brain, respectively. 

All the expression profiles of the mature miRNAs were 
normalized, and statistically significant differences 
between the profiles were identified as indicated above 
for the mRNA-seq and as previously described (20). For 
the comparison between the expression levels in the pineal 
and in the brain, a minimum fold change of 2 was required 
in addition to the above. 



Microarray and data analysis 

Labeled RNA preparation and hybridization to DNA 
microarrays were performed according to the Affymetrix 
manual with the two-cycle target labefing protocol (http:// 
www.affymetrix.com/support/downloads/manuals/expres 
sion_analysis_technical_manual.pdf). A total of 12 
Affymetrix DNA microarrays were hybridized with 
RNA-pools of pineal glands from fight and dark treat- 
ments in six time points throughout the daily cycle. 
Each DNA microarray was normalized using Affymetrix 
GeneGhip Operating Software. The entire DNA micro- 
array data set, logarithmically transformed, was 
normalized using quantile normalization to guarantee 
that the distribution of probe intensities was the same in 
aU the chips (24). 

The identification of transcripts with significantly differ- 
ent expression levels in the fight and dark groups was per- 
formed as follows: (i) The sum of differences in expression 
levels between the fight and dark groups was plotted against 
the transcripts average log expression for all transcripts rep- 
resented in the DNA microarray, along with a moving 
average trend line. Logarithmically transformed microarray 
data are highly variable at low expression levels, mainly due 
to the background (additive) noise (25). As a consequence, a 
moving average trend line is higher for low average expres- 
sion, (ii) Each transcript was given a z-score, measuring how 
far is its sum of differences from the trend fine, compared 
with the standard deviation of all the sum of differences of 
transcripts with a similar expression level. A high z-score for 
a given transcript indicates that the variation in the expres- 
sion levels is higher than the variation expected from the 
transcript average expression level. The specific z-score used 
as a detection cutoff was chosen to keep the associated false- 
positive rate at 5%. The false-positive rate for a given 
z-score was calculated by assuming a normal distribution 
for the transcripts expression levels, and accordingly, for 
the sum of differences in expression levels. Transcripts 
with low expression values (i.e. their mean expression level 
over all time points is lower than the average of all tran- 
scripts) in both the fight and dark groups were discarded. 
This analysis was performed for all six time points, yielding 
a fist of 16 genes. To avoid misdetection of genes due to 
inaccuracy in a single microarray measurement, the analysis 
was also performed for every five consecutive time points, 
yielding 13 additional genes (Supplementary Table SI). 

Gene Ontology analysis 

The 92 pineal gland light-induced genes that were detected 
by mRNA-seq (Table 1) were subjected for gene enrich- 
ment analysis using the DAVID functional annotation 
tool (26). The DAVID's default zebrafish genes back- 
ground was used. Significantly enricfied ontology terms 
(Benjamini-Hocfiberg false discovery rate of 5%) are 
fisted in Supplementary Table S2. 

Validation of light-induced gene expression by whole- 
mount in situ hybridization 

Wild-type zebrafish embryos were entrained by two I2-h: 
I2-h LD cycles and transferred to DD. On the third day of 
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Table 1. Pineal gland light-induced genes detected by mRNA-seq 



Number 


Gene syrnbol 


RefSeq ID 


L/D-fold 


Number 


Gene symbol 


RefSeq ID 


L/D-fold 


J 


nfl'fhd] 


iNiVl l700iU 


5 09 


47 


I'tiidi']/ fi 

1 l.tlHi lAU. 




2 28 


2 






4.92 


48 


V SI1 11 a 


NM 001007366 

i>iVl L/UiUU/JUU 


2.26 


3 




lS IVl UU 11/ LSji-y 


4 81 


49 


J/. L tiz 1 1-1 yjul J . 1 


MM 001077'^09 

IN iVl UU 1 U / / JUZ 


9 9S 

Z.Z J 


4 


^ VCllt} / 


IN IVl IX) J J I'-r 


4 81 


so 

JU 




MM lQ08fSQ 

IN iVl 1 770U7 


9 91 

Z.Z J 


5 


nrldl ( vcvcvha) 


INIVI Z,U J / Z7 


A 7Q 

ly 


51 


hi fin] 
II 1/ 1 til 


MM 90040S 

IN iVl ZUUH-U J 


2 22 


5 




NTM 00]00?Sfi'^ 

IN IVl UU 1 UUZ JO J 


4 61 


52 


-///Z JO 


MM 001089QS9 

IN iVl UU 1 U0Z7 JZ 


2 21 


7 


taghi2 


IN IVl IX) 1. J l\j 


4 SI 


53 


sfk35 


MM 001077f^lr^ 

IN iVl UU 1 U / / U 1 U 


2 2 


g 


nppc 


XTM 001 1 f^l "^41 

IN IVl UU 1 1 U 1 JH- 1 


4 38 


54 


eng2a 


MM 1 1 1 044 
IN ivi 1 J 1 yj'-t'-x 


2 2 


Q 


ffiinp 


iNivi iJioyj 


4 If, 

4-. ZD 


J J 


oactinl 


MM 1 1 0"^ 1 
IN iVl 1 J lU J 1 


9 1 7 
z. 1 / 




~gc 122979 


MM 0010^7^74 

IN IVl UU 1 U J / J / H- 


4 09 


56 


if I'd 1 


MM 00107^^SSS 

IN iVl UU 1 U / 0 J J J 


2 12 


1 1 

1 i 


j K op J 


MM 71 ^ 1 ilQ 
IN IVl Z 1 J 1 '4-7 


J.77 


S7 

J / 


hivep2 


MM 001 0'^OI f^A 

iNiVl UUIUJUIO'4- 


9 1 
z. 1 


12 




MM 100QS7 

IN IVl 1 yyyo 1 


0. 1 1 


58 


piKJlJU 


MM 901 14"^ 

IN iVl ZU 1 1 ^ J 


9 OQ 

Z.U7 


1 ^ 


Ufli/icfU f licC 1 J 


MM 717f^7Q 

IN IVl Z 1 ZU / 7 


3 6 


59 




MM 90SS^S'^ 

IN iVl ZU J JO J 


2 05 


14 


lutli J i-l 


MM 0017779'^4 

IN IVl UU 1 Z / / Z jH- 


-J. J7 


r^o 

UU 


IISlll 1 uz. 


MM 919790 

IN iVl Z 1 Z / ZU 


2 05 




oshplV 


MM 00100^097 
INIVI UUIUUJVZ/ 


J. j4 


(^1 

D 1 


11/7/? ^ fodhnd j 


MM 001 1 Q70S8 
iNiVl UUII7/UJ0 


9 0'^ 
Z.UJ 


1 f\ 




MM 91^'^zlQ 

IN IVl Z 1 J J'-ry 


^ 41 
j.t 1 


62 




MM 1898/^6 

IN iVl 1 oZoUO 


9 09 
z.uz 


1 7 

1 / 


njuj-o 1 CHopH—o J 


MM 001007918 
IN IVl UUlUUZZlo 


J.JO 


D J 


zgaojoj/ 


IN iVl Z 1 J J D J 


9 09 
Z.UZ 


1 R 
i o 


klf9 


MM 001 1 9879Q 

IN IVl UU 1 1 Zo / Z7 


J. J J 


64 


devind4a 


MM 001080Q8Q 

IN iVl UU 1 U0U7O 7 


9 01 

Z.U 1 


19 


itj'}i2ch 


MM 100Q80 

IN IVl 1 7 7 70U 


J. JZ 


65 


pptfilfib 


MM 001 1 O'^ 1 1 7 

IN iVl UU 1 1 U J 1 1 / 


2 




gcihvh2 


MM 001094'^R7 

IN IVl UU 1 UZH- Jo / 


^ 9Q 

J.Z7 


66 


coch 


MM 00100^89^ 

IN iVl UU 1 UU JoZ J 


1 99 


21 




MM OOIO'^OIR'^ 

IN IVl UU 1 U JU 1 0 J 


"K 9Q 

J.Z7 


67 


camK4 


MM 001017(^07 

IN iVl UU 1 U 1 / UU / 


1 96 


22 




MM 91^17Q 

IN IVl Z 1 J 1 / 7 


'K 98 

J.ZO 


68 


canih.2n 1 


MM 001 14S08Q 

IN iVl UU 1 1H-JU07 


1 ^fi 
1 .70 


23 


tuft la 


MM 001080001 

IN IVl UU 1 UoUUU 1 


J. 1 


69 


s}iQp25a. 


MM 1'^14'^S 

IN iVl 1 J 1 H- J J 


1 OS 
1 .7 J 


24 


pel 


MM 1898S7 

IN IVl 1 oZO J / 


'X 07 
J.U / 


70 


midn 


MM 9070S9 

IN iVl ZU / U JZ 


1 OS 

1 .7 J 




cry la 


MM 0010779Q7 

IN IVl UU 1 U / / Z7 / 




71 


c lac 


MM 001 1 1 01 9r^ 

IN iVl UU 1 1 1 u 1 ZU 


1 04 


26 


tillSTill VS 


MM 00107r^fS40 

IN IVl UU 1 U / UOH-U 


9 99 


72 


ucp 


MM 1 1 1 7^^ 

IN iVl 1 J 1 1 / 0 


1 01 
1 .7 1 


27 


tl/tiJ-Z f cHUpH- z, j 


MM 001 1 070^^S 

IN IVl UU 117/ UO J 


2 97 


73 


SCHo CIQ 


MM 1"^ 1(^98 

IN iVl 1 J 1 OZo 


1 91 


78 

Zo 


cci}iik2d2 


MM 001 009 '^49 
iNiVl UUIUUZJM-Z 


9 8Q 

Z.07 


74 
/ 4- 


c\'p27cl 


MM 001 1 1 '^'^'^7 
IN iVl UU 1 1 1 J J J / 


1 0 
1 .7 


29 


tlijpiJ { (AupA J 


MM 001 1 070(^9 

IN IVl UU 117/ UOZ 


9 87 
z.o / 


75 




MM 001080091 

IN iVl UU 1 UoUUZ 1 


1 78 
1 . / 0 


30 


sp tlc2h 


MM 001 1 14741 

IN IVl UU 1 1 1 H- / 4 1 


9 7*5 

Z. / J 


76 


myl9b 


MM 91^919 

IN iVl Z 1 J Z 1 Z 


1 lA 
1 . /*+ 


31 


pimS 


MM 0010'^4Q78 

IN IVl UU l\j j'-ry i o 


9 71 
z. / 1 


77 


rpc()5(i 


MM 9007S1 

IN iVl ZUU / J 1 


1 I'X 
1 . / J 


32 




MM 1 898/^9 

INIVI 1 CiZ.o\jZ. 


9 7 
z. / 


78 


sic 17 a? 


MM 0010087SS 

IN iVl UU 1 U70 / J J 


1 68 
1 .00 


"^"^ 


fund 


MM 001 1 9S'^49 
INIVI UUllZoJ'+Z 


9 'n7 
Z. J / 


7Q 

/ 7 


toil 


MM 1 'X 1 400 

IN iVl 1 J IM-UU 


1 f>l 
1.0/ 






MM 0010'^Q107 
IN IVl UUIUJ7IU/ 


9 '^4 
Z. J4- 


so 

oU 


grk 7ci 


MM OOIO'^I 841 
iNiVl UUlUJlo'4-1 


1 6 
1 .0 


35 


111 I J n J 


MM 7n4^Q 


2.54 


81 


till tfZ 1/ \ ' ^ Vc 1 UU ^ J 


1 > iVl 1 J 1 UU J 


1.6 


DU 


n pas 4 a 


MM 00104^^71 
INIVI UUIU'+JJZI 


9 S9 
Z. JZ 


89 
oZ 


ncLirod 


MM 1 "^0078 
IN iVl 1 JU7 / 0 


1 S8 
1 .Jo 


37 


S.fap 


NM 131373 


2.51 


83 


cplx4a 


NM 001077300 


1.55 


38 


speg 


NM 001007109 


2.5 


84 


slmn2a 


NM 001005923 


1.54 


39 


cry 3 


NM 131786 


2.48 


85 


alplalh 


NM 131690 


1.53 


40 


nrld2a (reverhhl) 


NM 001130592 


2.45 


86 


atplhla 


NM 131668 


1.45 


41 


ahcfl 


NM 201315 


2.4 


87 


dirasl 


NM 199831 


1.44 


42 


zgc:154093 


NM 001077731 


2.36 


88 


stxhpl 


NM 001025182 


1.43 


43 


argl 


NM 199611 


2.34 


89 


atp2h2 


NM 001123238 


1.41 


44 


clqtnf4 


NM 001017702 


2.33 


90 


gpm6aa 


NM 213200 


1.4 


45 


camk2hl 


NM 001271393 


2.29 


91 


gpm6ab 


NM 214687 


1.37 


46 


zgc:153615 


NM 001077791 


2.28 


92 


~ic2a 


NM 131558 


1.28 



development, embryos were exposed to a 1-h light pulse 
('daylight' fluorescent lamp, Ught intensity of 5W/m^) at 
six time points (CT2, 6, 10, 14, 18 and 22); control 
embryos remained in darkness. The collected embryos 
were fixed and subjected to whole mount in situ hybrid- 
ization (ISH) followed by quantification as previously 
described (11,27). DNA fragments from the coding 
regions of decl, e4bp4-6, perla, reverbhl, pfkfh4l and 
lonrfll were PCR amplified, cloned into pGEM T-Easy 
vectors (Promega) and served as templates for generation 
of anti-sense riboprobes using DIG RNA labeling kit 
(Roche). PCR amplification was conducted using the 
primers listed in Supplementary Table S3A. For each tran- 
script, vahdation of hght-induction by whole-mount ISH 
was conducted at two selected circadian time points. 
Differences between treatments were determined by two- 
tailed /-test. 



Validation of light-induced gene expression by 
quantitative real-time PCR 

Adult Tg(aa«a/2:EGFP) zebrafish were exposed to light 
and dark treatments as described in the section 'mRNA- 
seq, niiR-seq and microarray experimental procedures' 
above, resulting with three pools of six pineal glands per 
treatment under six sampling times. Total RNA was ex- 
tracted using RNeasy Lipid Tissue Mini Kit (Qiagen), and 
cDNA was prepared using High-Capacity cDNA Reverse 
Transcription Kit (Applied Biosystems). Quantitative 
real-time PCR (qRT-PCR) analysis was carried out by ab- 
solute quantification using ABI PRISM 7300 Sequence 
Detection System (Applied Biosystems), according to the 
manufacturer's protocol. Triplicates of light and dark 
cDNA samples from selected time points served as tem- 
plates for PCR amplification using SYBR Green PCR 
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Master Mix (Applied Biosystems) and specific primers 
(Supplementary Table S3B). Expression levels were 
normalized to beta-actin2 {bactin2) transcript, and differ- 
ences between treatments were determined by two-tailed 
?-test. Light-induction of every transcript was assayed by 
qRT-PCR at two circadian time points. 

Transient transfections of Pac-2 cells and live-cell 
bioluminescence assays 

Pac-2 cells (28) were cultured as described (29) and tran- 
siently transfected using the FuGene HD reagent (Roche 
Diagnostics) according to the manufacturer's protocol. 
Promoter reporter constructs: perlb-Luc contains 3.3 kb 
of perlb promoter region cloned into luciferase expression 
vector pGL3Basic [Promega, (30)]; cryla-Luc consists of a 
fragment of 1.3 kb encompassing 1.25 kb of the 5' flanking 
and 53 bp of cry la first exon, cloned into pGL3Basic (12); 
E-boXpcrih-Luc contains four copies of the perlb 
E-box (y-GAAG CACGTG TACTCG-30 cloned into the 
minimal promoter luciferase reporter vector pLucMCS 
[Stratagene, (30)]; D-boXp^,,.7-Luc contains six copies of 
the per2 D-box (5'-TGCGt "cTTATGTAAA AAGAG-3') 
cloned into pLucMCS (13); and AP-l-Luc reporter con- 
struct, consisting of four copies of the canonical crvla AP- 
1 #1 site (5'-TGACTCA-3') cloned into pLucMCS (12). 
The coding sequences of selected transcription factors 
were incorporated into CMV promoter-driven expression 
vectors: decl and reverbbl were cloned into pcDNA3.1 
(Invitrogen); e4bp4-5 and e4bp4-6 were cloned into 
pCS2-MTK [derivative of pCS2, (31)]. Empty expression 
vectors served as negative controls. Transiently trans- 
fected cells were exposed to three LD cycles and two 
24-h periods under DD conditions, and hve-ceU biolumin- 
escence assays were carried out as described elsewhere 
(13,30) using a Topcount NXT counter (2-detector 
model, Perkin Elmer). Traces represent the mean of four 
independently transfected wells ± SD from a single ex- 
periment. Each experiment was performed a minimum 
of two times. Peak and amplitude values (Supplementary 
Table S4) were calculated as in previous studies (13,32). 
Briefly, short-term and long-term trends were removed 
from the raw data by an adjacent-averaging method 
with 4-h and 5-day running means, respectively, providing 
'luciferase activity values'. Normalized values were 
obtained by dividing the luciferase activity values by 
their average value. The peak was calculated as the 
highest point of the normalized values on the third LD 
cycle. The amplitude was computed as the difference 
between the normalized luciferase activity of the peak 
and the preceding trough, divided by 2. Statistical differ- 
ences were determined by two-tailed /-test. 

Morpholino-mediated knockdown in embryos 

Knockdown experiments were performed using 
morpholino-modifled antisense oligonucleotides (MO, 
Gene Tools): Gene Tools standard control MO (5'- 
ctcttacctcagttacaatttata-3'); per2(E4i4)MO (5'-tgcagatgtact- 
tacagtgtttttg-3') designed to target the splice site between 
exon 4 and intron 4, resulting in aberrant per2 splicing 
(33); decl(E2i2)MO (5'-aattacgcattcgtaccttactgtc-3') 



targeting the exon 2 and intron 2 boundary to interfere 
with decl splicing; and miR-183(guide)MO (5'-acagtgaattc- 
taccagtgccatac-3'), designed to target dre-miR-183 guide 
strand and block its activity. One-cell stage embryos were 
injected with MOs (2 nl, 1 mM) and incubated under dif- 
ferent fighting regimens as indicated. The efficiency of 
per2(E4i4)MO and decl(E2i2)MO targeting splice sites was 
evaluated by RT-PCR ensuring the formation of altered 
transcripts throughout the locomotor activity monitoring 
period [5-8 days post-fertilization (dpf)]. Total RNA was 
extracted from injected embryos, and cDNA was 
synthesized as previously described (33). Transcripts 
were PCR-amplified using the primers listed in 
Supplementary Table S3C. 

Larval locomotor activity assays 

Embryos injected with control MO, per2(E4i4)MO 
or decl(E2i2)MO, were transferred within 4h post- 
fertilization (hpf) to a temperature-controlled incubator 
at 28° C and maintained under constant darkness for 
four days. On the fifth day, larvae were placed in 48-well 
plates in the observation chamber of the DanioVision 
Tracking System (Noldus Information Technology), en- 
trained by 3-h fight (LED, fight intensity of 1.8 W/m^) 
and transferred to constant dim light conditions 
(LED, light intensity of 0.013 W/m^) for four daily 
cycles. Locomotor activity was tracked and analyzed by 
the Ethovision 8.0 software (Noldus Information 
Technology). Tracks from days 6-8 were analyzed for 
the total distance moved (cm) by each larva per lOniin 
time-bins. The data are presented as a moving average 
(20 sfiding points) of ~24 larvae in each group. To deter- 
mine alterations in rhythmic locomotor activity, individ- 
ual tracks underwent Fourier analysis as previously 
described (34,35). Briefly, the time-dependent signals 
were converted to frequency-dependent signals using 
the Fast Fourier Transform and scored with a ratio 
('G-factor') representing the extent to which the frequency 
corresponds to a 24-h period. Differences in the G-factor 
distributions between control MO- and per2(E4i4)MO- or 
decl(E2i2)MO-treated groups were determined by the 
Kolmogorov-Smirnov test. 

Profiling light-induction of mRNA in Pac-2 cells by 
quantitative real-time PCR 

Pac-2 cells were exposed to three 12-h: 12-h LD cycles. 
Subsequently, total RNA was extracted from cells kept 
in darkness or exposed to fight [tungsten light source, 
0.2 W/m^, as previously described (12)] for different time 
periods (0.5, 1, 1.5, 2, 3, 4, 6 and 8h). Total RNA was 
isolated using Trizol reagent (Gibco BRL), and cDNA 
was transcribed using Superscript III reverse transcriptase 
(Invitrogen) with random primers. qRT-PCR was con- 
ducted using TaqMan Gene Expression Assays (Appfied 
Biosystems). PCR reactions were carried out using the 
ABI PRISM 7300 Sequence Detection System (Appfied 
Biosystems) with TaqMan Universal PCR Master Mix 
(Applied Biosystems) according to the manufacturer's 
instructions. Expression levels were calculated by the 
2-AACT jjjgjjjQ^ relative to bactin2. Values are the mean 
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of three independent cell pools. Differences between photic 
treatments and sampling times were determined by two- 
way analysis of variance (ANOVA). TaqMan pre- 
designed primer and probe sets (Assays-on-Demand, 
AppUed Biosystems) were used for decl (assay ID: 
Dr03101060_gl) and per2 (assay ID: Dr03093705_ml). 
Custom TaqMan (Assays-by-Design, AppHed 
Biosystems) was used to determine the expression of 
e4bp4-5, e4hp4-6, cryla and bactin2 with primers and 
probe sequences as presented in Supplementary 
Table S3D. 

Transient transfections and dual iuciferase assays in 
HEK-293T cells 

The constructs used for transient transfections of HEK- 
293T cells were prepared as follows. psiCHECK-aanat2 
and psiCHECK-e4bp4-6: the 3'UTRs of aanat2 and 
e4bp4-6 genes, containing binding recognition sites for 
niiR-183, were cloned downstream to the Renilla 
Iuciferase open reading frame in psiCHECK-2 Vector 
(Promega) that includes a control Firefly Iuciferase 
reporter under a different promoter. A 680 bp fragment 
of aanat2 3'UTR was PCR-amplifled from zebrafish 
embryo cDNA using specific primers (Supplementary 
Table S3E), both containing NotI restriction sites. The 
PCR product was digested with NotI and ligated into 
Notl-digested psiCHECK-2. A 2547 bp fragment of 
e4bp4-6 3'UTR was cloned into psiCHECK-2 m the 
same manner using specific primers (Supplementary 
Table S3E). psiCHECK-aanat2-M183 and psiCHECK- 
e4bp4-6-M183: site-directed mutagenesis was used to spe- 
cifically mutate 2nt of the miR-183 seed site in aanat2 
3'UTR and 3nt of the miR-183 seed site in e4bp4-6 
3'UTR, using QuikChange Site-directed mutagenesis kit 
(Stratagene) as instructed by the manufacturer. The sets 
of complementary primers used for generation of the 
desired point mutations are presented in Supplementary 
Table S3F. miRVec-183: the pre-miR-183 cloned into the 
niiRNA expression vector (miRVec) was provided by 
Prof. R. Agami (36). 

One day before transfection, HEK-293T cells were 
seeded in 96-well plates (Greiner Bio-One) at a density 
of 7500 cells per well, in Dulbecco's Modified Eagle 
Medium (Gibco BRL) supplemented with 10% fetal calf 
serum and 1% penicillin-streptomycin. Cells were trans- 
fected in quadruplicate using the jetPRIME Reagent 
(Polyplus-Transfection) according to the manufacturer's 
protocol, with 40 ng DNA: 4ng psiCHECK-2 enclosing 
wild-type or mutated 3'UTR or an empty vector; 36 ng 
miRVec containing pre-miR-183 or an empty vector. 
Control wells were transfected with 40 ng pEGFP-1 
vector (Clontech). 48 h after transfection, transfection ef- 
ficiency was assessed by GFP fluorescence in control wells, 
and Firefly and Renilla Iuciferase activities were measured 
using the Dual-Luciferase Reporter Assay System kit 
(Promega) and a Veritas microplate luminometer 
(Turner BioSystems), according to the manufacturer's in- 
structions. Co-transfection experiments were performed in 
three independent repeats. 



aanat2 mRNA rhythm assay 

The following procedure was designed to examine the 
effect of miR-183 knockdown on the aanat2 mRNA 
rhythm in the pineal gland. Embryos were injected with 
control MO or miR-183(aiiide)MO, exposed to two 12-h: 
12-h LD cycles and sampled with 4-h intervals under DD 
during days 3^ (CT2, 6, 10, 14, 18, 22 and 2b). Whole- 
mount ISH and quantification was conducted as formerly 
described (11) using the aanat2 probe (37). This procedure 
was repeated twice. 

RESULTS 

Characterization of the light-responsive coding 
transcriptome in the zebrafish pineal gland 

With the aim of identifying potential key players in light- 
entrainment of the circadian clock, we studied the effect of 
hght exposure on gene expression in the zebrafish central 
clock organ, the pineal gland. Light-induced phase shifts 
of the clock are known to be time-of-day-dependent (38); 
therefore, it is Hkely that clock-related hght-responses vary 
as a function of the circadian time. Taking this into 
account, the general experimental design included 1-h 
hght treatments at six time points, evenly distributed 
throughout the daily cycle (Figure 1 and 'Materials and 
Methods' section). Characterization of the light-induced 
pineal gland transcriptome was carried out both 
by mRNA-seq and microarray analyses. mRNA-seq 
analysis was performed on triplicates of hght and dark 
pineal gland samples that were pooled from six time 
points. Microarray analysis was conducted on six light 
and six dark pineal gland samples that were collected at 
4-h intervals, thus providing information regarding the 
hght-induction of genes at different circadian times. 

The mRNA-seq analysis yielded a hst of 92 light- 
induced genes in the pineal gland, which represent 
RefSeq-annotated genes (Table 1 and 'Materials and 
Methods' section). The microarray analysis revealed 29 
hght-induced RefSeq genes (Supplementary Table SI 
and 'Materials and Methods' section). The higher 
number of light-induced genes detected by mRNA-seq 
reflects the limited gene representation on the microarrays 
and the enhanced precision of mRNA-seq in measurement 
of transcripts levels (16). Importantly, 19 (65%) of 29 
RefSeq genes that were identified as light-induced by 
microarray were also detected by mRNA-seq. 

A subset of genes was selected for vahdation by 
qRT-PCR on adult pineal glands or by whole-mount 
ISH in embryos (Figure 2), repeating the light and dark 
treatments as in the original experimental procedure 
(Figure 1 and 'Materials and Methods' section). The ex- 
perimental validation of 9 of the 10 tested genes, together 
with the overlap between the two independent transcrip- 
tome methods, provides a high level of confidence for the 
identified hght-induced genes. 

Functional classification of light-regulated genes 

The global impact of light exposure on biological 
processes was analyzed by DAVID annotation enrichment 
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Figure 2. Validation of light-induction by whole-mount ISH and qRT-PCR. Upregulation of nine transcripts following 1-h of light (L) versus 
control darkness (D) is confirmed by quantitative whole-mount ISH on 48-72-h-old embryos (A-F) or by qRT-PCR on adult pineal gland samples 
(G-I). For ISH analysis, dorsal views of representative specimens are shown; signals in the pineal gland (A-C, green arrow) or the entire cranium 
area (D-F) were determined using ImageJ software. Error bars represent SE (A-F, n = 15-18; G-I, n = 3); *P<0.05, **P <0.0l, ***P <0.00\, 
t-test. CT, circadian time. 



tool (26). Significantly overrepresented Gene Ontology 
terms within the list of light-induced genes are listed in 
Supplementary Table S2 (Benjamini-Hochberg false dis- 
covery rate of 5%). Notably, multiple annotation terins 
associated with transcriptional regulation are significantly 
enriched (Supplementary Table S2), confirming the con- 
tribution of transcriptional control to the relay of photic 



information. An interesting finding is that ontology terms 
related to negative transcriptional regulation are highly 
represented (Supplementary Table S2), specifying the 
role of transcriptional repression in fight signafing. 

The enriched 'basic-leucine zipper domain' category 
supports a preferential role for basic-leucine zipper 
(bZlP) transcription factors in mediating light-regulated 
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gene expression. These include members of the PAR/ 
E4BP4 bZIP subfamily and AP-1 transcription factor 
complex that have been functionally linked to the 
process of hght-entrainment of the circadian clock 
(12,13,39). 

Interestingly, Gene Ontology terms related to the 
response to reactive oxygen species (ROS), such as 
'response to hydrogen peroxide', are significantly 
enriched, coinciding with the suggested involvement of 
ROS in hght signaling (40). 

Importantly, the most significantly enriched term is 'cir- 
cadian rhythm - mammal' that consists of seven genes that 
are considered as orthologs of mammalian clock genes, 
some of which were not previously recognized as light 
regulated (see below and Table 2). 

Circadian clock-related genes are induced by light 

Light-induced transcription of the core clock genes per2 
and cryla is well documented in most zebrafish tissues and 
cell lines and considered to play a fundamental role in the 
mechanism of clock entrainment by hght (10,11). As men- 
tioned above, additional clock gene orthologs that are 
induced by light in the pineal gland were identified by 
DAVID analysis. Furthermore, we identified several 
other hght-induced clock accessory loop genes, bringing 
the total number of hght-induced core clock and clock 
accessory loop genes in the pineal gland to 14 (Table 2). 
Besides per2 and cryla, two additional members of the 
negative limb of the core clock loop were encountered: 
cry3, which is also hght-induced in zebrafish Pac-2 cells 
(15) and per la. Transcriptional regulators belonging to 
clock accessory loops that were detected: decl and dec2, 
basic helix-loop-helix (bHLH) transcription factors that 
repress CLOCK: BMAL-mediated transcription (41); 
reverha, reverbbl and reverbb2, members of the nuclear 
receptor (NR) family that contribute to the robustness 
of the clock mechanism by regulating the expression of 
core clock genes (42); tefl, dbp2, e4bp4-2, e4bp4-5 and 
e4bp4-6 belonging to the PAR/E4BP4 bZIP subfamily of 
clock-controlled genes that feedback onto the core clock 
(13,43-45). Confirming our results is the fact that some of 
the pineal hght-induced genes (Table 2) were also 
identified previously, by microarray analyses, in other 
tissues or in cell hues (Supplementary Table S5). 

Cellular metabolic pathways are induced by light 

Hypoxia-inducible factor 1 alpha (hifla) is a key player in 
the regulation of cellular adaptations to hypoxia. It is a 
PAS-domain bHLH transcription factor that binds to core 
DNA sequence within the hypoxia response element of 
target gene promoters (46). Previous studies estabhshed 
a cross-talk between hypoxic signahng and circadian 
pathways. In one case, the expression levels of the core 
clock components PERI and CLOCK were increased in 
response to hypoxia in the mouse brain, an effect that was 
suggested to be mediated by direct interaction with HIFla 
(47). In another example, a hnk between hypoxic and cir- 
cadian pathways has been shown to exist in the combined 
transcriptional regulation of vasopressin by HIFla and 
CLOCK via an E-box enhancer (48). Recently, evidence 



Table 2. Circadian clock-related light-induced genes in the pineal 
gland, detected by mRNA-seq 
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analysis. 



for the interaction between the two pathways have been 
demonstrated in zebrafish, including direct binding 
of HIFla to E-boxes in the perl promoter, and both 
dampened perl circadian oscillations and loss of larval 
circadian activity patterns in response to hypoxia (49). 
In the present study, hiflal was identified and vahdated 
as being light-induced in the pineal gland (Table 1 and 
Figure 2H). 

Interestingly, a well-characterized target of HIFla (50), 
pfkjb4 {6-phosphofructo-2-kinase/fructose-2,6-biphosphatase), 
an isoenzyme involved in regulation of glucose metabolism 
in response to hypoxia, is highly regulated by light in the 
pineal gland and cranial areas (Table 1 and Figure 2E). This 
gene exhibits the greatest L/D-fold expression according to 
both transcriptome analyses. The finding that hiflal and its 
target are upregulated by light, together with the established 
cross-talk between hypoxic and circadian pathways, suggests 
that the hypoxic pathway is activated by photic signahng 
and contributes to circadian clock entrainment. 

Differential effects of light-induced transcription factors 
on clock-controlled and light-driven promoter activation 

The clock-related light-induced transcription factors, decl, 
reverbbl, e4bp4-5 and e4bp4-6, were selected for further 
analysis. The effects of these factors on clock- and light- 
regulated promoter activity were tested using the light- 
entrainable zebrafish cell line Pac-2 and a live-cell 
biolumiiiescence monitoring system. The perlb-Luc and 
E-boXp^^i-ifj-Luc constructs are markers for clock-controlled 
expression (30), whereas expression of the luciferase 
reporter constructs cryla-Luc and D-boXp„.2-Luc is pri- 
marily light-driven but with evidence for co-regulation 
by the circadian clock (12,13). Pac-2 cells were co- 
transfected with expression vectors for each of the tested 
transcription factors, together with the promoter/ 
enhancer-reporter constructs, exposed to different 
lighting regimens, and luciferase activity was monitored 
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and compared with that of control-treated cells (Figure 3 
and 'Materials and Methods' section). 

Based on the luciferase activity values, it is evident that 
over-expression of decl leads to a robust attenuation of 
the expression driven by the perlh promoter and E-box 
elements (Figure 3A). However, examination of the 
normalized luciferase activity levels reveals that the 
rhythmic expression patterns themselves were not abol- 
ished (Supplementary Figure SI). The attenuation of 
luciferase expression levels by decl was not observed for 
D-box elements or a control non-rhythmic AP-l-Luc 
reporter, and only a moderate reduction of expression 
was evident for the cryla promoter (Supplementary 
Figure S2). These effects of decl could be attributed to 
repression of CLOCK: BMAL-induced transcription via 
E-box enhancers (41). Furthermore, decl over-expression 
induced a phase delay of the rhythmic expression driven 
by the cryla promoter and D-box elements under LD 
cycles (P< 0.0001, /-test), such that the anticipatory rise 
and dechne of promoter activity were eliminated, as 
demonstrated by a peak at ~ZT12 (Figure 3 A and 
Supplementary Table S4). This effect is presumably due 
to the repression of the clock machinery by decl, leading 
to a predominantly light-induced expression pattern of the 
cryla promoter and D-box elements, with an increase in 
expression during the light phase and a decrease during 
the dark phase, thus lacking the clock-generated anticipa- 
tion of the hght-dark and dark-light transitions. 

Over-expression of reverbhl resulted in a decreased 
amplitude of the rhythms driven by the perlb promoter 
(/■< 0.0001, ?-test) and E-box elements (P<0.01, ?-test; 
Figure 3B and Supplementary Table S4), most likely due 
to repression of endogenous BMAL transcription by 
reverbhl through specific response element [RREs, (51)]. 
Moreover, reverbhl led to a phase shift of the rhythmic 
expression of D-box elements (/"< 0.001, /-test) and, to 
a lesser extent, also of cryla promoter-driven expression 
(P<0.01, /-test; Figure 3B and Supplementary Table S4). 
This outcome resembles the effect induced by decl, both 
resulting in a predominantly hght-driven regulation as a 
consequence of an inhibited clock. 

Interestingly, over-expression of e4bp4-5 and e4bp4-6 re- 
sulted in dissimilar effects. High levels of e4bp4-6 led to 
disruption of the expression driven by D-box elements 
(/■< 0.0001, /-test; Figure 3D and Supplementary 
Table S4), most hkely due to its direct binding to 
D-box enhancers and transcriptional repression (52). 
However, e4bp4-5 induced a lower amphtude of perlb 
promoter-driven rhythmic expression {P < 0.05, /-test; 
Supplementary Table S4), yet the expression pattern 
regulated by D-box elements remained rhythmic 
(Figure 3C). These results point to distinct functions for 
the two hght-induced e4bp4 homologs, in which e4hp4-6 is 
more related to light-regulation and e4hp4-5 to general 
clock function. 

Thus, constitutive over-expression of the four light- 
induced transcription factors resulted in differential 
effects on clock- and light-regulated expression in Pac-2 
cells, demonstrating various mechanisms by which tran- 
scription factors can mediate light-input to the core clock 
machinery and modulate clock gene expression. 



Light-induced genes are important for light-entrainment of 
rhythmic locomotor activity 

per2, which was identified by our transcriptome analyses, 
has been shown to play a key role in the onset of the pineal 
clock in vivo, demonstrated by its importance in the estab- 
lishment of the clock-controlled aanat2 mRNA rhythm 
that is triggered by light pulses (11,33). We further 
studied the role of per2 in vivo by testing the effect of 
per2 knockdown on the light-induced generation of 
clock-regulated behavior in larvae. Zebrafish larvae 
display circadian rhythms of locomotor activity 
characterized by a peak of activity during the subjective 
daytime (53). To entrain rhythms of locomotor activity, 
5 dpf larvae were exposed to 3 h of hght, and their loco- 
motor activity was subsequently monitored under 
constant dim light conditions. The timing of the light 
pulse (beginning of the day or night) determined the 
peak time of locomotor activity rhythms (day or night, 
respectively), indicating that the hght manipulation is re- 
sponsible for the initiation of clock-controlled rhythmic 
activity and its entrainment (Supplementary Figure S3). 
per2 knockdown [using per2 morphohno, per2(E4i4)MO] 
resulted in significantly disrupted patterns of circadian 
locomotor activity, as compared with control morpholino 
(cont MO)-injected embryos (P< 0.001; Figure 4A and C 
and 'Materials and Methods' section), demonstrating its 
role in hght-entrainment of clock-regulated behavior. 

We have shown that decl is induced by light in the adult 
and larval pineal gland (Table 1 and Figure 2A). Light- 
induction of decl has also been demonstrated in zebrafish 
larvae by microarray (15) and in the mouse SCN (41). The 
above experimental procedure was repeated to test the 
role of decl in the light-induced onset of clock-driven 
locomotor activity. Knockdown of decl [using decl 
morphohno, decl(E2i2)MO] significantly affected the circa- 
dian locomotor activity patterns triggered by hght pulses 
(P< 0.0001; Figure 4B and D and 'Materials and 
Methods' section). These results support our hypothesis 
that hght-induction of decl is important to the process 
of light-entrainment of the circadian clock. Together, 
these experiments demonstrate how the negative clock 
components, decl and per2, which act at the central 
clock organ by different mechanisms, contribute to the 
hght-entrainment of the circadian behavior. 

Kinetics of light-dependent changes in mRNA levels 

The kinetics of hght-induced expression might provide im- 
portant clues as to the position and function of genes in 
the hght signahng cascades. Therefore, we determined the 
detailed time course of light-induction of per2, cryla, 
e4hp4-5, e4bp4-6 and decl transcripts in the hght-respon- 
sive Pac-2 cells. Following entrainment by three LD 
cycles, Pac-2 cells were kept in darkness or exposed to 
light and sampled at eight time periods. Transcript levels 
were determined by qRT-PCR (Figure 5). The expression 
of per2, cryla and e4bp4-6 was significantly affected by 
light treatment (P< 0.0001, two-way ANOVA) and 
samphng time (P<0.01, two-way ANOVA). These tran- 
scripts display transient light-regulated expression 
patterns, in which mRNA levels rise following light 
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Figure 3. Light-induced transcription factors differentially affect clock- and light-regulated rhythmic promoter activity. Bioluminescence assay of 
zebrafish photosensitive Pac-2 cells transiently co-transfected with decl (A), reverbbi (B), e4bp4-5 (C) and e4bp4-6 (D), together with perlb-Luc (A- 
D, top left), E-boXp„]h-Luc (A-D, top right), cryla-Luc (A-D, bottom left) and D-hoXp„2-l^itc (A-D, bottom right). Normalized luciferase activity or 
luciferase activity values (counts per second, CPS) are plotted on the y-axis and time (days) is plotted on the x-axis. White and black bars above each 
panel represent the light and dark periods, respectively. Error bars represent SD. 
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Figure 4. Knockdown of perl and decl disrupts light-entrainment of rhythmic locomotor activity. Locomotor activity of 6-8 dpf larvae injected with 
control morpholino (cont MO; A and B, blue trendline), per2 morpholino [per2(i74i4)MO; A, red trendline] or decl morpholino [decl(E2i2)MO; B, 
orange trendline] under constant conditions. Injected larvae were raised in constant darkness and entrained by 3-h light pulse (indicated by arrows) 
on day 5, which determined the phase and peak time of rhythmic locomotor activity as displayed by control larvae. The average distance moved (cm/ 
lOmin) is plotted on the y-axis and zeitgeber time (ZT) is plotted on the x-axis. Error bars stand for SE (« = 24). Horizontal boxes indicate the 
lighting conditions before and throughout activity monitoring: Black boxes represent dark, white boxes represents light and diagonally lined boxes 
represent dim light. Significant differences in the distribution of the G-factors (a representation of the extent to which the frequency of the activity 
pattern corresponds to a 24-h period, see 'larval locomotor activity assays' in 'Materials and Methods' section) were revealed between the control 
morpholino- (C, black diamonds) and per2 morpholino (C, gray diamonds)-treated groups, as well as between the control morpholino- (D, black 
diamonds) and decl morpholino (D, gray diamonds)-treated groups {P < 0.001 and < 0.0001, respectively). The median G-factor value for each group 
is indicated (C-D, red lines). 
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Figure 5. Time course of light-induction in Pac-2 cells. Cells were exposed to three LD cycles. Subsequently, total RNA was extracted from cells 
kept in darkness (filled boxes) or exposed to light (empty boxes) for different time periods (0.5, 1, 1.5, 2, 3, 4, 6 and 8h). Quantification of per2 (A), 
cryla (B), e4hp4-6 (C), e4hp4-5 (D) and decl (E) niRNA levels was performed by qRT-PCR. Values shown are the mean of three independent cell 
pools; error bars represent SE. Statistical differences between photic treatments and sampling times were determined by two-way ANOVA. 
ZT = zeitgeber time, ZTO corresponds to the time of light onset. 



onset, reach maximal levels within a few hours and grad- 
ually descend, comparing with low basal levels under dark 
conditions (Figure 5A-C). e4bp4-6 mRNA is upregulated 
1.5 h after hght onset, shghtly before perl induction, 
peaks at 3h and gradually decays up to 6h after hght 
onset, similar to the per2 transcript (Figure 5A and C). 
Together with the finding that e4hp4-6 regulates the per2 
D-box (Figure 3D), these expression patterns suggest that 
e4bp4-6 regulates perl hght-induction. In contrast to 
e4bp4-6, the expression pattern of its homolog, e4hp4-5, 
is not light-regulated in Pac-2 cells (Figure 5D). 
Furthermore, unlike the case in the pineal gland, the 
decl transcript is also not induced by light in Pac-2 cells 
(Figure 5E), consistent with a previous report in the 
zebrafish BRF41 ceU line (54). Notably, the transient 



nature of the light-induction of genes, as demonstrated 
for per2, cryla and e4bp4-6, is most likely essential for 
their function. These transient patterns of light-induction 
imply the existence of regulatory mechanisms governing 
the kinetics of mRNA degradation. 

Pineal-enhanced and light-induced microRNAs 

One regulatory mechanism that may contribute to the 
transient nature of the hght-induced increase of specific 
transcripts involves targeting by miRNA. In the search 
for miRNAs that might be involved in circadian regula- 
tion and hght-entrainment, miR-seq was used to profile 
the miRNA population whose expression is enriched in 
the pineal gland and is hght-regulated ('Materials and 
Methods' section). Pineal glands were collected from 
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Table 3. The miR-183/96/182 cluster is pineal-enhanced and hght-induced 



miRNA 


Pineal D" 


Pineal L,^ 


L/D-fold 


Brain" 


Pineal/brain-fold 


dre-miR-182 


4339677 


7242855 


1.67 


2008 


2884 


dre-miR-96 


4108 


4735 


1.15 


5 


908 


dre-miR-183 


114631 


169606 


1.48 


178 


800 



"Normalized number of reads. 



light- and dark-treated fish at six circadian times 
(Figure 1) and pooled together, generating 'light' and 
'dark' samples. In addition, brains were collected from 
light- and dark-treated fish at CT2 and CT14 and 
pooled to generate a reference brain sample. 

This analysis resulted in the identification of 32 
niiRNAs whose expression is enhanced in the pineal 
gland relative to the brain (Supplementary Table S6), 
and 31 miRNAs that are upregulated by light in the 
pineal gland (Supplementary Table S7). The overlap of 
these two sets consists of 14 miRNAs that exhibit pineal 
enhanced and hght-induced expression (Supplementary 
Tables S6 and S7). Of particular interest is the miR-183/ 
96/182 cluster that is highly enriched in the pineal gland 
and is also upregulated by light (Table 3). These findings 
are consistent with the abundance of this cluster in the rat 
pineal gland and retina (55) and with previous reports of 
light-regulation of this cluster in the mouse retina (56), 
thus pointing to a conserved function in vertebrates. 

miR-183 targets the 3'UTR of the light-induced 
gene e4hp4-6 

To identify potential targets of the pineal-enhanced 
and hght-induced miR-183/96/182 cluster, we used 
TargetScanFish, a tool based on recognition site 
matching and features of site context [release 6.2, http:// 
www.targetscan.org; (57)]. Focusing on predicted target 
genes that are characterized by a transient light-induced 
expression pattern, we selected e4bp4-6, which encloses a 
putative binding site for miR-183 in its 3'UTR, for further 
analysis. To test the ability of miR-183 to downregulate 
e4bp4-6 mRNA levels, we used the luciferase reporter 
assay ('Materials and Methods' section). The 3'UTR of 
e4bp4-6, either with a wild-type sequence or carrying mu- 
tations in the miR-183 binding site (Figure 6B), was 
cloned downstream to a Renilla luciferase reporter gene. 
Reporter plasmids and miR-183 expression vector were 
co-transfected into HEK-293T ceUs. The relative activity 
of the Renilla luciferase reporters was normalized to 
control Firefly luciferase reporters and compared with 
control transfections with an empty miRNA expression 
vector. miR-183 significantly reduced the Renilla 
\uc\iemse-e4bp4-6 activity to 49%, compared with the 
control vector. However, when the miR-183 binding site 
was mutated the effect was statistically insignificant 
(Figure 6A). These results suggest that miR-183 may 
directly downregulate e4bp4-6 mRNA levels by interacting 
with a target site in its 3'UTR. Furthermore, this finding 
serves as an example for the possible role of miRNAs in 




■ miRVec 
□ miR-183 



no 3'UTR e4bp4-6 3'UTR e4bp4-6 mut 



B 

e4bp4-6 3'UTR 1 4 8 3-GC7\AGAAAGAUCAAGGUGCCAUA... 

1 1 1 1 1 1 1 

miR-183 3' -GUCACUUAAGAUGGUCACGGUAU-5' 

I I 1 1 

e4bp4-6 mut 1483-GCAAGAAAGAUCAAGGAGGCACA... 

Figure 6. miR-183 regulates the e4hp4-6 3'UTR in HEK-293T cells. 
(A) miR-183 represses the luciferase activity of a Renilla reporter 
that harbors the wild-type e4hp4-6 3'UTR but not with mutated 
miR-183 target site {e4hp4-6 mut) or empty reporter vector (no 
3'UTR). Luciferase activity is normalized to the activity of Firefly 
luciferase co-expressed from the dual reporter and compared with a 
co-transfection of a negative control miRNA vector (miRVec). A rep- 
resentative experiment is shown. Error bars represent SE; ***/*< 0.001, 
r-test with Bonferroni correction. (B) Representation of the miRNA- 
mRNA interaction tested in the luciferase assay, including the mutated 
miR-183 binding site. 



the generation of transient induction of transcripts in 
response to light, by promoting mRNA degradation. 

niiR-183 targets the 3'UTR of the clock-controlled gene 
aanat2 

Subsequently, we asked whether light-induced 
miRNAs could also contribute to the generation of 
mRNA rhythms of clock-controlled genes. aanat2 
(arylalkylamii2e-N-acet}'ltransferase-2), the key enzyme in 
melatonin synthesis (58), is specifically expressed in the 
zebrafish pineal gland in a clock-dependent manner (37), 
and contains a potential target site for miR-183. We tested 
the direct interaction of miR-183 and its putative target 
site in aanat2 3'UTR using the luciferase reporter assay, as 
described in the previous section. A significant reduction 
of the Renilla luciferase activity (57%, compared with the 
control vector) was generated by miR-183 when the 
aanat2 3'UTR was included, an effect that was abohshed 
when the miR-183 target site was mutated (Figure 7). 
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Figure 7. miR-183 regulates the aanat2 3'UTR in HEK-293T cells. 
(A) miR-183 represses the luciferase activity of a RenUla reporter that 
harbors the wild-type aanat2 3'UTR but not with mutated miR-183 
target site (aanat2 mut) or empty reporter vector (no 3'UTR). 
Luciferase activity is normalized to the activity of Firefly luciferase 
co-expressed from the dual reporter and compared with a 
co-transfection of a negative control miRNA vector (miRVec). 
A representative experiment is shown. Error bars represent SE; 
***/>< 0.001, /-test with Bonferroni correction. (B) Representation of 
the miRNA-mRNA interaction tested in the luciferase assay, including 
the mutated miR-183 binding site. 



This indicates that miR-183 might directly control aanat2 
mRNA levels through binding a target site in its 3'UTR, 
and adds support to the proposed involvement 
of miRNAs in the regulation of circadian mRNA 
rhythms (59). 

miR-183 regulates pineal aanat2 mRNA rhythms in vivo 

To test whether miR-183 regulates aanat2 in vivo, we 
studied the effect of miR-183 morphohno-mediated 
knockdown on the rhythms of pineal gland aanat2 
mRNA levels. Zebrafish pineal aanat2 transcription 
displays a robust circadian rhythm starting at 2 dpf (37). 
Morpholino-injected embryos were maintained under two 
LD cycles to entrain the circadian clock and then 
transferred to constant darkness. During the third and 
fourth days of development, embryos were collected at 
4-h intervals and subjected to whole-mount ISH for 
aanat2 mRNA. aanat2 mRNA levels were significantly 
affected by morpholino treatment and sampling time 
(P< 0.001, two-way ANOVA). Specifically, miR-183 
knockdown resulted in elevated aanat2 mRNA levels 
compared with the control treatment at the end of the 
subjective night and beginning of subjective day (CT22 
and CT2), time points at which a dechne of aanat2 
mRNA normally occurs (Figure 8). These results suggest 
a role for miR-183 in the regulation of aanat2 mRNA 
stability in vivo and its contribution to the rhythms of 
aanat2 mRNA in the pineal gland. 




Figure 8. miR-183 knockdown effects the pineal aaiuitl mRNA 
rhythm. Zebrafish embryos were injected with control morpholino 
(cont MO) or miR-183 morphohno [miR-183(gj,i£ie)MO] and entrained 
by two LD cycles. On the third day of development, embryos were 
transferred to DD, sampled at 4-h intervals and pineal aanatl 
mRNA levels were determined by whole-mount ISH. Gray and black 
horizontal bars represent subjective day and night, respectively. Values 
represent the mean ± SE optical densities of the pineal signals 
(n = 15-20). Representative pineal aanat2 signals (dorsal view) of 
embryos treated with cont MO (top panel) and miR-183(,,uije)MO 
(bottom panel) from each group are presented. **P<0.01, 
***/*< 0.001, /-test with Bonferroni correction. CT, circadian time. 



DISCUSSION 

The zebrafish constitutes a fascinating vertebrate model 
for studying the regulation of the circadian clockwork 
by light (6). With the goal of further exploring the mech- 
anisms underlying light-entrainment of the vertebrate 
circadian clock, we profiled the light-induced coding tran- 
scriptome and miRNome in the zebrafish central clock 
organ, the pineal gland, using cutting edge genome-wide 
methods. 

An interesting outcome of this approach is the identifi- 
cation of 14 core clock and clock accessory loop genes that 
are induced by light in the pineal gland (Table 2). The 
finding of light-regulation of multiple clock components 
in a central clock tissue points to a complex regulatory 
mechanism underlying hght-entrainment of the circadian 
clock. Notably, with the exception of two genes (tefl and 
dhp2), all these light-induced circadian clock-related genes 
function as transcriptional repressors, suggesting that 
light-entrainment is primarily mediated through the 
negative limbs of the oscillator. A similar result using 
microarray technology has been reported in the light- 
induced transcriptome of zebrafish larvae, heart and cell 
culture, in which six of the eight light-regulated clock 
genes that were identified act as transcriptional repressors 
(15). The functional analysis of selected light-regulated 
accessory loop genes {decl, reverhbl, e4bp4-5 and e4bp4- 
6) yielded differential effects on clock- and light-regulated 
expression in Pac-2 ceUs (Figure 3), demonstrating various 
mechanisms by which transcription factors modulate 
clock gene expression and thereby convey photic informa- 
tion to the core clock machinery. 

The circadian clock regulator, decl, strongly attenuated 
the clock-controlled expression driven by the perlb 
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promoter and E-box elements without affecting tlie phase 
(Figure 3A), probably reflecting its repression of 
CLOCK: BMAL-induced transcription (41). This repres- 
sion of the clock machinery by decl possibly underlies 
the elimination of the clock-derived anticipation in the 
cryla- and D-box-driven expression, resulting in purely 
hght-induced expression patterns (Figure 3A). We have 
shown that decl is highly induced by light in the adult 
and larvae pineal gland (Table 1 and Figure 2A). This is 
also the case in the mouse SCN (41), but not in zebrafish 
cell lines [Figure 5E, (54)]. A study in rfec7-deficient mice 
provided evidence for its role in circadian clock resetting 
(60). Our data indicate a similar role for decl in the 
process of light-entrainment in zebrafish; knockdown of 
decl resulted in disrupted circadian locomotor activity 
patterns that were triggered by a single light pulse 
(Figure 4B and D), resembling the effect generated by 
per2 knockdown (Figure 4A and C). Light-induction of 
per 2 and decl is predominant in the pineal gland. 
Nevertheless, whether the effects of these genes on entrain- 
nient of rhythmic behavior reflect their transcriptional 
effects and activities in the pineal gland or in peripheral 
clocks remains to be elucidated. 

The D-box C7>acting element has been shown to serve 
a central role in fight-driven gene expression in 
zebrafish (12,13,15). We have previously identified and 
characterized the zebrafish PAR/E4BP4 D-box binding 
factors (27). Here we show that two e4bp4 homologs, 
e4bp4-5 and e4hp4-6, both are induced by light in the 
pineal gland (Table 1 and Figure 2B and I), but only 
e4bp4-6 is induced by light in Pac-2 cells (Figure 5C and 
D). Furthermore, these factors yielded differential effects 
on clock- and Hght-regulated promoter activation in Pac-2 
cells (Figure 3C and D). In particular, e4hp4-6 but not 
e4bp4-5, abofished expression driven by D-box elements, 
presumably by direct binding and transcriptional repres- 
sion (52). Further studies will be required to assess the 
relative contribution of the various D-box-binding bZlP 
factors to the fight response. 

The transient nature of light-induced gene expression 
implies to the involvement of mechanisms that control 
mRNA stabifity, such as regulation by miRNA. In the 
search for candidate miRNAs that may play a role in 
light-entrainment and circadian regulation, the pineal 
gland-enhanced and fight-induced miRNAs were profiled 
by miR-seq (Supplementary Tables S6 and S7). This 
analysis identified the miR-183/96/182 cluster, which is 
highly enriched in the pineal gland and upregulated by 
light (Table 3). This miRNA cluster has been formerly 
shown to exhibit a diurnal variation of expression in the 
mouse retina, and evidence were provided for their target- 
ing of adcy6, a clock-controlled gene that modulates 
melatonin synthesis (61). In another study in the mouse 
retina, the miR-183/96/182 cluster was found to be hght- 
regulated and to target the voltage-dependent glutamate 
transporter slclal (56). This cluster has also been shown 
to be abundantly expressed in the rat retina and pineal 
gland, where it displays daily dynamics of expression 
(55). Together, these findings suggest a conserved 
function for the miR-183/96/182 cluster in vertebrates 
and support its involvement in circadian and light- 



regulation. We have demonstrated that miR-183 
downregulates e4bp4-6 mRNA through a specific target 
site in its 3'UTR (Figure 6). Moreover, we have shown 
that miR-183 targets the 3'UTR (Figure 7) and regulates 
the rhythmic mRNA profile of the pineal clock-controlled 
gene, aanat2, in vivo (Figure 8). These findings emphasize 
the contribution of miRNAs to the transient expression 
profiles of light- and clock-regulated genes, and their 
importance to pineal functions, such as fine-tuning the 
kinetics of rhythmic melatonin production. Moreover, 
we have generated a database of additional miRNAs 
that may serve as candidates for regulating features of 
the vertebrate circadian clock and its entrainment by light. 

The identity of the photopigments that mediate light- 
induction in the zebrafish pineal gland remains unclear, 
although several possible candidates have been suggested 
[as recently reviewed (6)]. Exo-rhodopsin as well as several 
other extra retinal opsins are expressed in the pineal gland 
and may act as photoreceptors (62-64). Another candi- 
date opsin photoreceptor is teleost multiple tissue 
(TMT) opsin, whose mRNA is detected in the central 
nervous system, most peripheral tissues and even in cefi 
fines (65) and has been confirmed to function as a photo- 
receptor (66). Consistently, most zebrafish tissues and 
even cell lines are directly fight-responsive (8). CRYs 
were also suggested to act as photopigments due to the 
fact that CRY plays a role in photoreception by lateral 
neurons in Drosophila (67), and given that one of the 
zebrafish cry homologs {cry4) shares higher sequence simi- 
larity with Drosophila cry than with the mammalian 
homologs (68). Finally, fight has been shown to induce 
the production of ROS such as H2O2 that might fink 
photoreception to the induction of the zebrafish core 
clock genes cryla and per2 (40). It is suggested that 
fight-activated enzymes such as flavin-containing oxidase 
function as photoreceptors that trigger accumulation of 
intracellular H2O2 upon exposure to near violet-blue 
wavelengths of light, as demonstrated in zebrafish Z3 
cells (40). In addition, an interesting finding of the 
current study is that hifla, as well as its characterized 
target pjkfb4l, are induced by light (Table 1 and Figure 
2E and H). Together with observations of molecular 
cross-talk between hypoxic signaling and circadian 
pathways (47^9,69), these findings hint at the existence 
of an additional photic signaling cascade that might be 
involved in the process of light-entrainment. Any of the 
above photic transduction pathways may contribute to 
activating a set of fight-induced genes in the pineal 
gland. The characterization of the fight-induced transcrip- 
tome of the central clock organ points to the involvement 
of multiple pathways in the relay of photic information 
from photoreceptors to core circadian clock function. 

The approach taken in this study enabled the identifi- 
cation of light-induced genes that may be involved in 
fight-entrainment, the characterization of their action 
on the molecular clock and elucidating their effects on 
rhythmic behavior in the intact animal. Therefore, this 
study strengthens the utihty of the zebrafish as a model 
for exploring the basic mechanisms whereby circadian 
clocks respond to the environment and how they have 
adapted during evolution. 
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